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Abstract. The Laplace operator in infinite quantum waveguides (e.g., a bent strip 
or a twisted tube) often has a point-like eigenvalue below the essential spectrum that 
corresponds to a trapped eigenmode of finite L2 norm. We revisit this statement for 
resonators with long but finite branches that we call "finite waveguides". Although 
now there is no essential spectrum and all eigenfunctions have finite L2 norm, the 
trapping can be understood as an exponential decay of the eigenfunction inside the 
branches. We describe a general variational formalism for detecting trapped modes 
in such resonators. For finite waveguides with general cylindrical branches, we obtain 
a sufficient condition which determines the minimal length of branches for getting 
a trapped eigenmode. Varying the branch lengths may switch certain eigenmodes 
from non-trapped to trapped states. These concepts are illustrated for several typical 
waveguides (L-shape, bent strip, crossing of two stripes, etc.). We conclude that 
the well-established theory of trapping in infinite waveguides may be incomplete and 
require further development for being applied to microscopic quantum devices. 
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1. Introduction 

By its name, a waveguide serves for propagating waves which may be of different 
physical origins: fluctuations of pressure in acoustics, electromagnetic waves in optics, 
particle waves in quantum mechanics, surface water waves in hydrodynamics, etc. The 
transmission properties of a waveguide can be characterized by its resonance frequencies 
or, equivalently, by the spectrum of an operator which describes the waves motion (e.g., 
the Laplace operator in the most usual case). For infinite waveguides, the spectrum 
consists of two parts: (i) the essential (or continuous) spectrum for which the related 
resonances are extended over the whole domain and thus have infinite L 2 norms, and 
(ii) the discrete (or point-like) spectrum for which the related eigenfunctions have finite 
L 2 norms and thus necessarily "trapped" or "localized" in a region of the waveguide. 
A wave excited at the frequency of the trapped eigenmode remains in the localization 
region and does not propagate. 

The existence of trapped, bound or localized eigenmodes in classical and quantum 
waveguides has been thoroughly investigated (see reviews [H [2] and also references in 
[3]). In the seminal paper, Rellich proved the existence of a localized eigenfunction in 
a deformed infinite cylinder [I]. His results were significantly extended by Jones [5]. 
Ursell reported on the existence of trapped modes in surface water waves in channels 
0,0 IE], while Parker observed experimentally the trapped modes in locally perturbed 
acoustic waveguides [9j[T0]. Exner and Seba considered an infinite bent strip of smooth 
curvature and showed the existence of trapped modes by reducing the problem to 
Schrodinger operator in the straight strip, with the potential depending on the curvature 
[TT] . Goldstone and Jaffe gave the variational proof that the wave equation subject to 
Dirichlet boundary condition always has a localized eigenmode in an infinite tube of 
constant cross-section in any dimension, provided that the tube is not exactly straight 
[T2] . The problem of localization in acoustic waveguides with Neumann boundary 
condition has also been investigated [13], HI] . For instance, Evans et al. considered 
a straight strip with an inclusion of arbitrary (but symmetric) shape [H] (see [15] for 
further extension). Such an inclusion obstructed the propagation of waves and was 
shown to result in trapped modes. The effect of mixed Dirichlet, Neumann and Robin 
boundary conditions on the localization was also investigated (see 0, EE] and references 
therein). A mathematical analysis of guided water waves was developed in [17] . 

All the aforementioned works dealt with infinite waveguides for which the Laplace 
operator spectrum is essential, with possible inclusion of discrete points. Since these 
discrete points were responsible for trapped modes, the major question was whether or 
not such discrete points exist below the essential spectrum. It is worth noting that the 
localized modes have to decay relatively fast at infinity in order to guarantee the finite L 2 
norm. But the same question about the existence of rapidly decaying eigenfunctions may 
be formulated for bounded domains (resonators) with long branches that we call "finite 
waveguides" (Fig. [T|) . This problem is different in many aspects. Since all eigenfunctions 
have now finite L 2 norms, the definition of trapped or localized modes has to be revised. 
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Quite surprisingly, a rigorous definition of localization in bounded domains turns out to 
be a challenging problem [181 HH1 120] • In the context of the present paper concerning 
finite waveguides, an eigenmode is called trapped or localized if it decays exponentially 
fast in prominent subregions (branches) of the bounded domain. The exponential decay 
of an eigenfunction in the branch can be related to the smallness of the associated 
eigenvalue in comparison to the cut-off frequency, i.e. the first eigenvalue of the Laplace 
operator in the cross-section of that branch |21j . In other words, the existence of a 
trapped mode is related to "smallness" of the eigenvalue, in full analogy to infinite 
waveguides. Using the standard mathematical tools such as domain decomposition, 
explicit representation of solutions of the Helmholtz equation and variational principle, 
we aim at formalizing these ideas and providing a sufficient condition on the branch 
lengths for getting a trapped mode. The dependence of the localization character on the 
length of branches is the main result of the paper and a new feature of finite waveguides 
which was overseen in the well-established theory of infinite waveguides. As in practice 
all quantum waveguides are finite, this dependence may be important for microelectronic 
devices. 

The paper is organized as follows. In Sec. EJ we adapt the method by Bonnet- 
Ben Dhia and Joly [17] in order to reduce the original eigenvalue problem in the 
whole domain to the nonlinear eigenvalue problem in the domain without branches. 
Although the new problem is more sophisticated, its variational reformulation provides 
a general framework for proving the trapping (or localization) of eigenfunctions. We 
use it to derive the main result of the paper: a sufficient condition ( fl9|) on the 
branch lengths for getting a trapped mode. In sharp contrast to an infinite (non- 
straight) waveguide of a constant cross-section, for which the first eigenfunction is always 
trapped and exponentially decaying [12] , finite waveguides may or may not have such 
an eigenfunction, depending on the length of branches. This method is then illustrated 
in Sec. |3]for several finite waveguides (e.g., a bent strip and a cross of two stripes). For 
these examples, we estimate the minimal branch length which is sufficient for getting 
at least one localized mode. At the same time, we provide an example of a waveguide, 
for which there is no localization for any branch length. We also construct a family 
of finite waveguides for which the minimal branch length varies continuously. As a 
consequence, for a given (large enough) branch length, one can construct two almost 
identical resonators, one with and the other without localized mode. This observation 
may be used for developing quantum switching devices. 

2. Theoretical results 

For the sake of clarity, we focus on planar bounded domains with rectangular branches, 
while the extension to arbitrary domains in MJ 1 with general cylindrical branches is 
straightforward and provided at the end of this Section. 

We consider a planar bounded domain D composed of a basic domain Q of arbitrary 




Figure 1. Two examples of a finite waveguide: (a) a planar bounded domain D which 
is composed of a basic domain Q of arbitrary shape and three rectangular branches 
Qi of lengths and width 6=1; (b) a three-dimensional bounded domain with three 
general cylindrical branches. 

shape and M rectangular branches Qi of lengths a; and width b as shown on Fig. [TJ 

M 

D = ttu\jQ l . 

i=i 

We denote Tj = dfl D dQi the inner boundary between the basic domain Q and the 
branch Qi and T = dQ\ \J^ =1 Tj the exterior boundary of Q. We study the eigenvalue 
problem for the Laplace operator with Dirichlet boundary condition 

-AU= \U, U\ dD = 0. (1) 



2.1. Solution in rectangular branches 

Let Ui(x,y) denote the restriction of the solution U(x,y) of the eigenvalue problem (JT]) 
to the branch Qi. For convenience, we take 6 = 1 and assume that the coordinates x 
and y are chosen in such a way that Qi = {(x, y) e M 2 : < x < a i: < y < 1} (the 
final result will not depend on this particular coordinate system). The eigenfunction 
Ui(x, y) satisfying Dirichlet boundary condition on dQi has the standard representation: 

oo 

Ui(x, y) = U\ Q . = c n sinh(7 n (ai - x)) sm(irny), (2) 

71=1 



where 7„ = y ir 2 n 2 — A and c n are the Fourier coefficients of the function U at the inner 
boundary Tj (at x = 0): 

i 

J dy U(0,y)sin(7rny). (3) 



sinh(7 n a 

o 

Substituting this relation into Eq. (j5J) yields 

/ \ „^, 7r . , xx sinh(7„(aj — x)) . , . 

Ui(a;,3/) = 2 2^(^|r i ,sni(7m3/)) £a(ri) — s sm(7rm/), (4) 

sinn^7nQ-v 

where the integral in was interpreted as the scalar product in L 2 (Fj). The 
representation flj]) is of course formal because its coefficients are still unknown. 
Nevertheless, one can already distinguish two different cases. 
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(i) If A < 7r 2 , all 7„ are real, and the representation (j3J decays exponentially. In 
fact, writing the squared L 2 -norm of the function Ui(x, y) along the vertical cross-section 
of the branch Qi at x, 

1 



u 2 {x, y)dy = 2 ^(t7|r iS sin(7rny))^ (r . 

n=l 



sinh (n/njaj - x)) 
sinh 2 (7 n ai) 



one can use the inequality sinh(7„(aj — x)) < sinh^a^e llX to get 

oo 

/,(a;)<2^(f/| ri ,sin(7™y)) 2 2(ri)e - 2 ^ = / l (0)e- 2 ^ (0 < x < a t ). (5) 

n=l 

This shows the exponential decay along the branch with the decay rate 271 = 2 V ir 2 — X. 

(ii) In turn, if A > tt 2 , some 7 n are purely imaginary so that sinh (7^) becomes 
sin(|7 n |z), and the exponential decay is replaced by an oscillating behavior. 

One sees that the problem of localization of the eigenfunction in the basic domain 
Q is reduced to checking whether the eigenvalue A is smaller or greater than it 2 (or ir 2 /b 2 
in general). 

2.2. Nonlinear eigenvalue problem 

The explicit representation (pE} of the eigenfunction in the branch Qi allows one to 
reformulate the original eigenvalue problem ([1]) in the whole domain D as a specific 
eigenvalue problem in the basic domain Q. In fact, the restriction of U onto the basic 
domain Q, u = Up, satisfies the following equations 

-Au = Xu in ft, u\r = 0, u\ Ti = Ui\ Ti , 7^-|r\ = --^1^, (6) 

where d/dn denotes the normal derivative directed outwards the domain. The last 
two conditions ensure that the eigenfunction and its derivative are continuous at inner 
boundaries Tj (the sign minus accounting for opposite orientations of the normal 
derivatives on both sides of the inner boundary). The normal derivative of Ui can 
be explicitly written by using Eq. @: 

~dn^ T ' = ~~dx^ x= ° = 2 ^^ coth (lnai)(U lri ,sm(7rny)) L2{ri) sm(7Tny). (7) 

n=l 

Denoting Tj(A) an operator acting from if 1 / 2 ^) to H~ l l 2 (Ti) (see [22] for details) as 

00 

Ti(X)f = 2^7 n coth(7 n aj)(/, sin(7rm/)) L2(ri) sin(7rny), 



n=l 



the right-hand side of Eq. (JTj) can be written as 

The eigenvalue problem ()6]) admits thus a closed representation as 

Ou 

-Au — Xu mQ, u\r — 0, — | r , = -Ti(X)u\ ri . (8) 
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The presence of branches and their shapes are fully accounted for by the operators Tj 
which are somewhat analogous to Dirichlet-to-Neumann operators. 

Although this domain decomposition allows one to remove the branches and get 
a closed formulation for the basic domain Q, the new eigenvalue problem is nonlinear 
because the eigenvalue A appears also in the boundary condition through the operators 
Tj(A). A trick to overcome this difficulty goes back to the Birman-Schwinger method 
[231 [23] ( see a l so [2H]). Following [17], we fix A and solve the linear eigenvalue problem 

du 

- Au = fi(X)u inQ, u\ r = 0, —\ Ti = ~T i {\)u\ Tv (9) 

where //(A) denotes the eigenvalue which is parameterized by A. The solution of the 
original problem is recovered when /i(A) = A. 

From a practical point of view, a numerical solution of Eqs. (Q with the subsequent 
resolution of the equation /x(A) = A is in general much more difficult than solving the 
original eigenvalue problem (see also [26] for possible numerical schemes). In turn, Eqs. 
(J9]) are convenient for checking whether the first eigenvalue Ai is smaller or greater than 
7r 2 , as explained below. 

2.3. Variational formulation 

We search for a weak solution of the eigenvalue problem fl9]) in the Sobolev space 

Hl{ti) = {v(x,y) e L 2 (n), dv/dx e L 2 {Q), dv/dy E L 2 {Q), v\ r = 0}. 
Multiplying Eq. ([9]) by a trial function v G Hq(Q) and integrating by parts, one gets 

/i(A) J vu = — J vAu = J (Vf , Vm) — J v ~q~- 

n n n an 

Since v vanishes on T, the weak formulation of the problem reads as 

M 

(Vtz,V«) La ( n) +X;(^(^«)£ 2 (r 1 )=MA)(w,«) ia(n) VveH^Q). (10) 
The first eigenvalue /xi(A) is then obtained from the Rayleigh's principle 

M (v^v«) Mn , + S?f,(r.(A) iV ,) Mr ^ 

t»efl3(0), vjto ( v i v )L 2 (n) 
One can show that /ii(A) is a continuous monotonously decreasing function of A on 
the segment (0, 7r 2 ]. For this purpose, one first computes explicitly the derivative of the 
function 



h(X) = 7 n coth(7 n aj) = V n 2 n 2 — X coth(V ir 2 n 2 — X a^) 

and checks that h'(X) < 0. Now one can show that /ii(Ai) < /ii(A2) if Ai > A2. If some 
trial function v 2 minimizes the Rayleigh's quotient (fTTJ) for A2, one has 

/ii(A0 < 



(Vv a , Vw 2 ) L2(c) + Y.f=i( T i^i)v2,v 2 ) L2{Vi) 



{Vv 2 ,Vv 2 ) L2m +J2Zi(Ti{X 2 )v 2 ,v 2 ) L2ir . ) 
S 7 ^ — /^i(A 2 J, 

{V2,V 2 ) L2in) 
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where the monotonous decrease of the function h(X) was used (the mathematical proof 
of the continuity for an analogous functional is given in [27]). 

Since the function /xi(A) is positive, continuous and monotonously decreasing, the 
equation /ii(A) = A has a solution < A < 7r 2 if and only if fii(n 2 ) < it 2 . This is a 
necessary and sufficient condition for getting a trapped mode for the linear eigenvalue 
problem (jUJ). 

2.4- Sufficient condition 

For any trial function v G Hq(Q), we denote the Rayleigh's quotient 

fi(v) = L2(n) , ~ L2{Tl) . (12) 



Since 7 n (^ 2 ) = ir^/n 2 — 1, one has 7i(vr 2 ) = 0, and the operators Ti(ir 2 ) can be 
decomposed into two parts so that 

M 

i=i 1 

oo M 

+ 2n Vn 2 — 1 coth(7rat V n 2 — 1) (p , sin^ny))^^ j. (13) 

n=2 i=l 

If one finds a trial function p G #o(^) f° r which yu(p) < it 2 , then the first eigenvalue 
Hi{tt 2 ) necessarily satisfies this condition because ^i(tt 2 ) < ^(f)- The inequality //(p) < 
7r 2 is thus a sufficient (but not necessary) condition. Given that coth(7ra i v / ^ 2 — 1) < 
coth(7raj\/3) for any n > 2, the sufficient condition can be written as 

M AI 

J^— </3-^K i coth(7ra iV / 3), (14) 



a 

i=i 



where 



/3 =7r 2 (p,p) L2(n) -(Vp,Vp) L2(n) , (15) 



cii = 2(p,sin(vry)) i2(ri) , (16) 

oo 
n=2 

Before moving to examples, several remarks are in order. 

(i) The shape of the branches enters through the operator Tj(A). Although the above 
analysis was presented for rectangular branches, its extension to bounded domains in M. n 
with general cylindrical branches is straightforward and based on the variable separation 
(in directions parallel and perpendicular to the branch). In fact, the Fourier coefficients 
(u, sin(7rnp))L 2 (r i ) on the unit interval (i.e., the cross-section of the rectangular branch) 
have to be replaced by a spectral decomposition over the orthonormal eigenf unctions 
{■0n(y)K£Li of the Laplace operator Aj_ in the cross-section of the studied branch (in 
general, is a bounded domain in M™ -1 ): 

A±1pn + Vnlpn = m T h IpnldTi = 0. (18) 
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In particular, the operator Tj(A) becomes 

oo 

T iWf = ^lnCoth(-f n ai)(f,ilj n ) L2{ri) ilj n (y), 



n=l 



with 7„ = i/z/ n — A. Repeating the above analysis, one immediately deduces a sufficient 
condition for getting a trapped mode: 



M M 

a. 



^ — < - ^ /tiCoth(q i A/z/ 2 - z/i), (19) 



1=1 i=i 



with 



=Mv,v) L2(a) -(Vv,Vv) L2[a) , (20) 



= (vMl^), (21) 

oo 
n=2 

One retrieves the above results for rectangular branches by putting ijj n (y) = a/2 sin(7rra/) 
and z/ n = 7r 2 ?2 2 . 

The inequality (|T9l is the main result of the paper. Although there is no explicit 
recipe for choosing the trial function v (which determines the coefficients 0, and Ki), 
this is a general framework for studying the localization (or trapping) in domains with 
cylindrical branches. 

(ii) If the branches are long enough (e.g., a« ^> (u 2 — z/i) -1 / 2 ), the values 
coth(a iv /z/ 2 — ui) are very close to 1 and can be replaced by 1 + e where e is set by 
the expected minimal length so that the inequality ( JT9l) becomes more explicit in terms 
of a{. 

M M 

^</3-(l + e)5>. (23) 

8=1 1=1 

In the particular case when all <7j are the same, one can introduce the threshold value r\ 

as 

El /J (1 + e) >r-^ /nAS 

i=i «=i 
For domains with identical branches, = a, the above condition determines the branch 
length ath = M/ rj which is long enough for the emergence of localization. This means 
that for any a > a t h there is a localized eigenmode. Since a t h was obtained from the 
sufficient condition ( fT9l) . the opposite statement is not true: for a < a t h, this condition 
does not indicate whether the eigenfunction is localized or not. In fact, a th is the upper 
bound for the minimal branch length a m \ n which may distinguish waveguides with and 
without localized modes (see Sec. |3]). 

(iii) The trial function should be chosen to ensure the convergence of the series in 
( )22l) . If the boundary of Q is smooth, the series in (1221 converges for any function v from 
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Figure 2. Three types of a bent waveguide: (a) L-shape, (b) bent strip, and (c) 
truncated L-shapes parameterized by the length £ varying from (triangular basic 
domain) to 1 (the original L-shape). 



Hq(Q) according to the trace theorem [22] ■ In turn, the presence of corners or other 
singularities may require additional verifications for the convergence, as illustrated in 
Sec. E3 

(iv) The implementation of various widths 63 of the rectangular branches is relatively 
straightforward (e.g., sin(7rm/) is replaced by sm(7iny /bi), etc.). In order to guarantee 
the exponential decay in all branches, one needs A < tt 2 /&? for all i, i.e., A < 7r 2 / max{5 2 }. 
Rescaling the whole domain in such a way that max{6j} = 1, one can use the above 
conditions. 

(v) According to Eq. (jSJ), the decay rate, 2^, is determined by the eigenvalue A. 
Since jj,{y) is an upper bound for the first eigenvalue, one gets a lower bound for the 
decay rate: 

(M 
- ^[cxi/ai + m coth(aiy/u 2 - Vi)} 
i=i 

3. Several examples 

As we already mentionned, there is no general recipe for choosing a trial function v from 
Hq(Q). Of course, the best possible choice is the eigenfunction on which n(v) reaches 
its minimum. Except for few cases, the eigenfunction is not known but one can often 
guess how it behaves for a given basic domain. Since the gradient of the trial function v 
enters into the coefficient /3 with the sign minus, slowly varying functions are preferred. 
In what follows, we illustrate these concepts for several examples. 

3.1. L-shape 

We start by a classical problem of localization in L-shape with two rectangular branches 
of lengths a\ and a 2 (Fig. [2|i) for which the basic domain is simply the unit square. 
In the limit case a\ = a 2 = (i.e., D = Q, without branches), the first eigenvalue 
Ai = 2tt 2 > 7r 2 so that, according to our terminology, there is no localization. Since Ai 
continuously varies with a (01 = a 2 = a), the inequality A x > n 2 also remains true for 
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relatively short branches. In turn, given that Ai < it 2 for infinitely long branches, there 



should exist the minimal branch length a m ; n such that Ai 



7T 



At this length, the first 



eigenfunction passes from non- localized state (a < a m i n ) to localized state (a > a m i n ). 
In what follows, we employ the sufficient condition (fl~4l) in order to get the upper bound 
for a m \ n . 

The most intuitive choice for the trial function would be the first eigenfunction for 
the unit square with Dirichlet boundary condition, v(x,y) = sin(7rx) smijry). However, 
one easily checks that (3 = for this function, while a, and K{ are always non-negative. 
As a consequence, the condition (fT4|) is never satisfied for this trial function. It simply 
means that the first choice was wrong. 

For the trial function 



v (x, y) — (1 + x) sin(7n/) + (1 + y) sin(7r:r), 
the explicit integration yields 



(25) 



lv,v 



L 2 (n) 



dx I dy[(l + x) sin(7n/) + (1 + y) sin(7nc)]' 



-i 



-i 
o 



(Vu, Vv) 



dx / dy 



dx) \dy 



dv 



71 

T + 1 - 



1 2 

3 + 



(v,sm(-Ky)) L2{Ti) = 
(v,sm(imy)) L2(Ti} 



dy[(l + 0) sin(7T|/) + (1 



,sin(7rnx)) L2(r2) = 



y) sin(7r0 
(n > 2). 



sin(Try) = -, 



from which 







O"! = a 2 = -? Kl = K 2 = 0. 



The condition (TTJl) reads as 



1 1 

— + — < 2. 

ai a 2 



(26) 



If the branches have the same length, ai = a 2 = a, then the upper bound of the minimal 
branch length for getting a localized eigenfunction is given by a t h = 1. 

We also solved the original eigenvalue problem (CQ) for L-shape with a\ = a 2 = a by 
a finite element method (FEM) implemented in Matlab PDEtools. The first eigenvalue 
Ai as a function of the branch length a is shown by solid line on Fig. [3j One can 
clearly see a transition from non-localized (Ai > 7r 2 ) to localized (Ai < tt 2 ) states when 
a crosses the minimal branch length a m j n ~ 0.84. As expected, the theoretical upper 
bound a t h which was obtained from a sufficient condition, exceeds the numerical value 
a min . In order to improve the theoretical estimate, one has to search for trial functions 
which are closer to the true eigenfunction. At the same time, a t h and a m ; n are close to 
each other, and the accuracy of the theoretical result is judged as good. Similar results 



for L-shape in three dimensions are derived in Appendix A 
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Figure 3. The first eigenvalue Ai divided by it 2 , as a function of the branch length 
a (pi = a 2 = a), for three bent waveguides shown on Fig. [2 L-shape (solid line), bent 
strip (dashed line) and truncated L-shape with £ — (dash-dotted line). For the first 
two cases, the curves cross the level 1 at a m i n « 0.84 and a min « 2.44, respectively. 
In turn, the third curve always remains greater than 1 (see explanations in Sec. 13.41) . 
For a — 0, Ai is respectively equal to 2tt 2 , Aa 2 and 5tt 2 , where a ss 2.4048 is the first 
positive zero of the Bessel function Jq(z). 



3.2. Cross 

Another example is a crossing of two perpendicular rectangular branches (Fig. H^), 
for which the basic domain is again the unit square. Since the trial function ( )25|) also 
satifies the boundary condition for this problem, the previous sufficient condition ( 1261) 
remains applicable for arbitrary lengths 03 and 04. This is not surprising because any 
increase of the basic domain (i.e., if the basic domain was considered as the unit square 
with two branches Q3 and Q4) decreases the eigenvalue. A symmetry argument implies 
that other consecutive pairs of branch lengths can be used in the condition ( )26|) . e.g., 
the eigenfunction is localized if 1/02 + 1/03 < 2 for arbitrary a\ and a±. In turn, the 
condition 1/ai + 1/03 < 2 is not sufficient for localization (in fact, taking a 2 = a 4 = 
yields a rectangle without localization). 

The specific feature of the cross is that the exterior boundary of the basic domain 
Q consists of 4 corner points. We suggest another trial function 

v{x,y) = x(l + x)+y(l + y), (27) 

which satisfies the Dirichlet boundary condition at these points. The direct integration 
yields 

11 2 2 64 

^ 90 3' 7T 6 ' 



00 / l_(_l)n\ 2 

2tt V Vn 2 - 1 2 ^ 0.0029. 

^— ' V 7r 3 n 3 / 

r>=9 V / 



n=2 
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Figure 4. (a) Crossing of two rectangular branches; (b) an extension of the related 
basic domain f2; and (c) coupling between two waveguides from Fig. [2}: (£ = 0) with 
an opening of size e. 



The condition 



m 



reads now as 



E 1 



< 



ft 



«1 



coth(7raj\/3). 



(28) 

=1 ~* ' " " " 8=1 

If all the branches have the same length a, the upper bound of the minimal branch 
length can be estimated by solving the equation 
4 4ki 



0"! 



coth(7ra t h\/3), 



a th cti 

from which one gets ath ~ 0.2407. Note that this result proves and further extends 
the prediction of localized eigenmodes in the crossing of infinite rectangular stripes 
which was made by Schult et al. by numerical computation [28J. In that reference, the 
importance of localized electron eigenstates in four terminal junctions of quantum wires 
was discussed. 

Figure [5] presents first eigenfunctions for the crossing of two rectangular branches 
with aj = 5 (the second eigenfunction, which looks similar to the third one, is not shown). 
As predicted by the sufficient condition (|28|) . the first eigenfunction (with Ai ~ 0.667T 2 ) 
is clearly localized in the basic domain and exponentially decaying in the branches. In 
turn, all other eigenfunctions (with X n > ir 2 ) are not localized. 

It is worth noting again that any increase of the basic domain (see Fig. HJd) reduces 
the eigenvalue and thus favors the localization. 



3. 3. Bent strip 

In previous examples, the basic domain was the unit square. We consider another shape 
for which the analytical estimates can be significantly advanced. This is a sector of the 
unit disk which can be seen as a connector between two parts of a bent strip (Fig. [2b). 
In contrast to the case of infinite stripes for which Goldstone and Jaffe have proved the 
existence of a localized eigenmode for any bending (except the straigh strip) [12], there 
is a minimal branch length required for the existence of a localized eigenmode in a finite 
bent strip. In order to demonstrate this result, we consider the family of trial functions 

v a (r) = S -— (0 < a < 1). (29) 
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n = 1 n = 3 n = 4 n = 5 



Figure 5. First eigenfunctions for the crossing of two rectangular branches {a% = 5). 
The associated eigenvalues are Ai ss 0.6617T 2 , A2 = A3 ps 1.0327T 2 , A4 ps 1.036tt 2 and 
A 5 w 1.044tt 2 . 



In |Appendix B[ we derive Eqs. (IB.lt EH IB.3j) for the coefficients 0, ai, and «j 



respectively. Since all these coefficients depend on a, its variation can be used to 
maximize the threshold 77 given by Eq. ( 124)) . The numerical computation of these 
coefficients suggests that 77 is maximized for a around 1/3: 77 ~ 0.7154. If a% = 02 = a, 
one gets the upper bound a t h of the minimal branch length which ensures the emergence 
of the localized eigenfunction: 

2 

a > a t h = - ~ 2.7956. 
77 

We remind that this is sufficient, not necessary condition. The numerical computation 
of the first eigenvalue in the bent strip (by FEM implemented in Matlab PDEtools) 
yields a min 2.44. One can see that the upper bound a t h is relatively close to this 
value. The behavior of the eigenvalue Ai as a function of the branch length a is shown 
by dashed line on Fig. [3j 



3.4- Waveguide without localization 

Any increase of the basic domain Q reduces the eigenvalues and thus preserves the 
localization. In turn, a decrease of Q may suppress the trapped mode. For instance, the 
passage from L-shape (Q being the unit square) to the bent strip (Q being the quarter of 
the disk) led to larger minimal length required for keeping the localization (a min ps 2.44 
instead of a m i n ~ 0.84). For instance, when a\ = 02 = 2, the first eigenfunction, which 
was localized in the L-shape, is not localized in the bent strip (Fig. [6]). Further decrease 
of the basic domain Q may completely suppress the localization. 

In order to illustrate this point, we consider the truncated L-shape shown on Fig. 
|2b with I = for which 

n = {(x,y)eM. 2 : - 1 < x < 0, - 1 < y < 0, x + y > -1} 
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Figure 6. The first eigenfunction for three bent waveguides shown on Fig. [2](£ = 0), 
with a — 2. The associate eigenvalue Ai is equal to 0.93 5 77T 2 , 1.00867T 2 and 1.14357T 2 , 
respectively. Although the first eigenmode is localized, all three eigenfunctions visually 
look similar. 

is a triangle. It is easy to see that u(x, y) = cos(7ra;) + cos (717/ ) is the first eigenfunction 
of the following eigenvalue problem in Q: 

du 

—Au = hu in Q, u\r = 0, — |r- = 0, 

on 

with the eigenvalue fix = it 2 . From the variational principle 

(Vv, Vv) L2{n) 



fix = inf 



so that 



(Vv, Vv) L2{n) > fa(v,v) L2(n) = n 2 (v,v) L2(n) \/v G H%(0). 
Moreover, the Friedrichs-Poincare inequality in the branches Qi implies [22] 

(Vv, Vv) L2(Qi) > 7r 2 (v,v) L2m Wv G HfcQi), 

from which 

(Vv, Vv) L2{D) > n 2 (v, v) L2{D) Wv G Hq(D). 

As a consequence, all eigenvalues of the original eigenvalue problem (JT|) in D exceed 7r 2 
and the corresponding eigenfunctions are not localized in the basic domain Q, whatever 
the length of the branches. 

When one varies continuously I in Fig. |2b, the basic domain transforms from the 
unit square (Fig. |2k) to triangle so that one can get any prescibed minimal length a m i n 
between 0.84 and infinity. In other words, for any prescribed branch length, one can 
always design such a basic domain (such £) for which there is no localized eigenmodes. 
As a consequence, the localization may be very sensitive to the shape of the basic domain 
and to the length of branches. These effects which were overseen for infinite waveguides, 
may be important for microelectronic applications. 

Figure [7] shows the first eigenfunction for three bent waveguides from Fig. [2] with 
a = 20. The associate eigenvalue Ai is equal to 0.9302tt 2 , 0.9879tt 2 and 1.0032tt 2 , 
respectively. Although the last two values are very close to each other, the behavior 
of the associated eigenfunctions is completely different. According to the sufficient 
condition, the first two eigenfunctions are localized in the basic domain, while the last 
one is not. One can clearly distinguish these behaviors on Fig. [71 
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Figure 7. The first eigenfunction for three bent waveguides shown on Fig. 0(^ = 0), 
with a = 20. The associate eigenvalue Ai is equal to 0.9302vr 2 , 0.9879?r 2 and 1.00327T 2 , 
respectively. Although the last two values are very close to each other, the behavior 
of the eigenfunctions is completely different. 

3.5. Two coupled waveguides 

The coupling of infinite waveguides has been intensively investigated [29]. We consider 
a coupling of two finite crossing waveguides through an opening of variable size e as 
shown on Fig. [4b. When e = 0, one has two decoupled waveguides from Fig. [2b for 
which we proved in the previous subsection the absence of localized eigenmodes. When 
e = v2, there is no barrier and the waveguides are fully coupled. This is the case of 
crossing between two rectangular branches as shown on Fig. HH, for which we checked 
the existence of localized eigenmodes under weak conditions ( 1261) or ( l28i) . Varying the 
opening e from to a/2, one can continuously pass from one situation to the other. 
This transition is illustrated on Fig. [8] which presents the first eigenfunction for two 
coupled waveguides shown on Fig. [4b, with = 5 and different coupling (opening e). 
The first two eigenfunctions, with e = (fully separated waveguides) and e = 0.4a/2 
(opening 40%), are not localized, while the last two eigenfunctions, with e = 0.5-^/2 
(opening 50%) and e = y2 (fully coupled waveguides, i.e. a cross), are localized. The 
critical coupling e c , at which the transition occurs (i.e., for which Ai = it 2 ) lies between 
40% and 50%. Although numerical computation may allow one to estimate e c more 
accurately, we do not perform this analysis because the value e c anyway depends on the 
branch lengths. In general, for any a > a m j n ~ 0.84, there is a critical value e c (a) for 
which Ai = 7r 2 . For e < e c , there is no localized modes, while for e > e c there is at least 
one localized mode. The high sensitivity of the localization character to the opening e 
and to the branch lengths can potentially be employed in quantum switching devices 
(see [30J and references therein). 

Conclusion 

We investigated the problem of trapped or localized eigenmodes of the Laplace operator 
in resonators with long branches that we called "finite waveguides". In this context, 
the localization was understood as an exponential decay of an eigenfunction inside the 
branches. This behavior was related to the smallness of the associated eigenvalue A 
in comparison to the first eigenvalue of the Laplace operator in the cross-section of 



Trapped modes in finite quantum waveguides 

k- *- i 

e = (0%) e = 0.4^ (40%) £ = 0.5^2 (50%) 

Figure 8. The first eigenfunction for two coupled waveguides shown on Fig. 
with a* = 5 and different coupling (opening e): e = (fully separated waveguides, zero 
coupling), £ = 0.4V2 (opening 40% of the diagonal), £ = 0.5V2 (opening 50% of the 
diagonal) and e — y/2 (fully coupled waveguides, no barrier). The associate eigenvalue 
Ai is equal to 1.05tt 2 , 1.02tt 2 , 0.977T 2 , and 0.677T 2 , respectively. In the first two cases, 
the eigenmodes is not localized. Changing the opening e, one passes from non-localized 
to localized eigenmodes. 

the branch with Dirichlet boundary condition. Using the explicit representation of an 
eigenfunction in branches, we proposed a general variational formalism for checking the 
existence of localized modes. The main result of the paper is the sufficient condition 
( 1191) on the branch lengths for getting a trapped mode. In spite of the generality of 
the formalism, a practical use of the sufficient condition relies on an intuitive choice of 
the trial function in the basic domain (without branches). This function should be as 
close as possible to the (unknown) eigenfunction. Although there is no general recipe for 
choosing trial functions, one can often guess an appropriate choice, at least for relatively 
simple domains. 

These points were illustrated for several typical waveguides, including 2D and 3D 
L-shapes, crossing of the rectangular stripes, and bent stripes. For all these cases, 
the basic domain was simple enough to guess an appropriate trial function in order 
to derive an explicit sufficient condition for getting at least one localized mode. In 
particular, we obtained the upper bound of the minimal branch length which is sufficient 
for localization. We proved the existence of a trapped mode in finite L-shape, bent 
strip and cross of two stripes provided that their branches are long enough, with an 
accurate estimate on the required minimal length. These results were confirmed by a 
direct numerical resolution of the original eigenvalue problem by finite element method 
implemented in Matlab PDEtools. The presented method can be applied for studying 
the localization in many other waveguides, e.g. smooth bent strip [3], sharply bent strip 
[31] or Z- shapes [32] . 

It is worth emphasizing that the distinction between localized and non-localized 
modes is much sharper in infinite waveguides than in finite ones. Although by definition 
the localized eigenfunction in a finite waveguide decays exponentially, the decay rate 
may be arbitrarily small. If the branch is not long enough, the localized mode may be 
visually indistinguishable from a non-localized one, as illustrated on Fig. [61 In turn, the 
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Figure Al. L-shape in three dimensions for which the basic domain O = [—1, 0] 3 is 
the unit cube. 



distinction between localized and non-localized modes in infinite waveguides is always 
present, whatever the value of the decay rate. 

The main practical result is an explicit construction of two families of waveguides 
(truncated L-shapes on Fig. [2b and coupled waveguides on Fig. Hb), for which the 
minimal branch length a min for getting a trapped mode continuously depends on the 
parameter i or e of the basic domain. For any prescribed (long enough) branch length, 
one can thus construct two almost identical finite waveguides, one with and the other 
without a trapped mode. The high sensitivity of the localization character to the shape 
of the basic domain and to the length of branches may potentially be used for switching 
devices in microelectronics. 
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Appendix A. L-shape in three dimensions 

As we mentioned at the end of Sec. [21 an extension to other types of branches is 
straightforward. We illustrate this point by considering the L-shape in three dimensions, 
i.e. two connected parallelepipeds of cross-section in the form of the unit square, for 
which the basic domain Q is the unit cube (Fig. IA1I) . For each branch, the eigenvalues 
and eigenfunctions in Eq. ([18]) for the cross-section can be parameterized by two indexes 
m and n: v m ^ n = n 2 (m 2 + n 2 ) and tfj mjn (y,z) = 2 sm(7rmy) sin(7mz) (similar for the 
second branch). 

We take the trial function 



which satisfies the Dirichlet boundary condition. The coefficients (3, ai and can be 
found from Eqs. ( [201 12T1 1221) . for which the explicit integration yields 



v(x, y, z) = [(1 + x) sin(7T?/) + (1 + y) sin(7ra;)] sin(7T2;), 



o 












1 



1 



-1 
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/dv\ 2 /dv\ 2 /dv\ 2 ~\ 7r 2 3 
AdxJ \dyJ \dz, 



3 2' 



-l -l -l 



so that (3 = 2tt 2 (Vv } Vv] 



i a (n) - M)z 2 (n) = 1/2- 





^i-i^ 2 (rx) 



dy / u(0, y, z) 2 sin(7ry) sin(7rz) 



-l -l 





y dz t>(0,y,z) 2sin(7rmy) sin(Trnz) = 



-1 -1 

(for m ^ 1 or n ^ 1), from which <7i = 02 = 1/4 and K\ = K2 = 0. The condition (TT9]) 
reads as 

i + i<2. 

ai a2 

If the branches have the same length, a\ = a 2 = a, then the upper bound of the 
minimal branch length for getting a localized eigenfunction is given by a t h = 1, as in 
two dimensions. 



Appendix B. Computation for bent strip 



The computation of the coefficients and oi is straighforward, while that for Ki requires 
supplementary estimates. 
Coefficient (3. One has 

1 



71 



71 



— It "sin (Trr)dr = — 



4 



1 



2(1 -a) 



where 



w v (q) = \ r v cos(qr)dr. 



Similarly, 



(Vv 



7T 



r(v' a ) 2 dr 



71 



7i cos Txr a sin nr 



,1+a 



dr. 



Expanding the quadratic polynomial and integrating by parts, one gets 



(Vf„, Vv a ) L2 



(n) 



7i 

T 



I W 2 a-l(27T) 



Combining this term with the previous result yields 

7r 3 2a 



A 2a -I 



2(1 -a) 1 -2a 

Uf 

U7 2q -i(27t) 



(B.l) 
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In the limit a — > 1/2, one has 



lim W2 - l(2?r) = ^ « 0.2257, 

a-H/2 2a — 1 27T 

where Si(x) is the integral sine function. 

Coefficients crj. For these coefficients, one gets 

i 



v 



sin(7rr)) L2(ri) = J r a sm 2 (nr)dr = \ (yY~ ~~ w ^( 2n ) 



from which 



1/1 . ,y 

0"i = cr 2 = - w a {2n) 

I v 1 — a / 



Coefficients One considers 



(f , sin(7rnr)) 



i 2 (ri) 



o 
1 

2 L 



r a sin(7rr) sin(7rw)dr 



w a (n(n - 1)) - w a (n(n + 1)) 



The function w a (q) can be decomposed into two parts [33] . 

oo oo 

w a (q) = / cos(qr)dr — / r~ a cos(qr)dr = q 



2 a r(f) 



(B.2) 



- w a (g), 



where 



w a (q) = / r a cos(qr)dr. 



We have 



where 



2vr Vn 2 - 1 (v, Bin(7rnr)) = 2?r V^ 2 - 1 (4 - e„) 



n=2 



2 1+a r(f) 



n=2 



(n - l) a_1 - (n + 1) 



a-l 



e n = - w a (-n{n - 1)) - w a {-n(n + 1)) . 

In order to estimate the coefficients e n , the function w a (q) is integrated by parts 
that yields for q — ir(n ± 1): 

*a+2(g) 



COS 

iB a (g) = a — 5 a (a + 1)- 



oc- 



cosg 



cosg 



a(a + l)(a + 2) — -±- + a(a + l)(a + 2) (a + 3) 



U7 a+4 (g) 
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The inequality 

oo 

(a + 3)\w a+4 (q)\ < (a + 3) J r'^dr = 1 

i 

leads to 

cosg , lW cosg+1 _ cosq . cosg — 1 
a — 5 aaa + l)(o; + 2) < w a {q) < a — a{a + lHa + 2) . 

from which 

where the lower and upper bounds are 

± _ 1 r/ a a(a + l)(a + 2) (-l)™" 1 =F 1 
6n ~ 2L\^2(n- l) 2 tt 3 (n- l) 4 

/ a (-l) n+1 a(a + l)(a + 2) (-l) n+1 ± 1 
~ 1^2 ( n + 1)2 ^4 (n+ l) 4 

Using these estimates, one gets 

oo 
n=2 

oo oo 

^ V^ 2 - 1 rfne n > ^ Vn 2 - 1 d n e~ = Aj, 

n=2 n=2 

oo oo 

£ v^T e 2 < £ y/rf=l (e+) 2 ee 



n=2 n=2 



from which 



«i < re, « = 2tt(Ai - 2^2 (B.3) 

Although the expressions for Al, and A% are cumbersome, the convergence of these 
series can be easily checked, while their numerical evaluation is straightforward. 

The numerical computation of these coefficients shows that the threshold value rj 
is maximized at a « 1/3: rj rs 0.7154. Note that if the coefficients Kj were computed 
by direct numerical integration and summation, the value of 77 for a = 1/3 could be 
slightly improved to be 0.7256. The difference results from the estimates we used, and 
its smallness indicates that the estimates are quite accurate. 
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